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Abstract 

We discuss the advantages of combining the experimental bound on Br(B s — > 
and the measured Br(£> — > Kl + £~) to get the model independent constraints 
on physics beyond the Standard Model. Since the two decays give complementary 
information, one can study not only the absolute values of the Wilson coefficients 
that are zero in the Standard Model, but also their phases. To identify the sector in 
which the new physics might appear, information about the shapes of the transverse 
asymmetries in B — > K*£ + £~ at low q 2, s can be particularly useful. We also empha- 
size the importance of measuring the forward-backward asymmetry in B — > K£ + £~ 
decay at large q 2, s. 
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1 Introduction and basic formulas 



Ever since the first observation of B — > K*^ [1], the decay modes governed by the loop- 
induced b — > s transition play a major role in the search of signals of physics beyond 
the Standard Model (BSM) in low energy experiments. After years of experimental and 
theoretical effort we now know that the observed decay rates and several conveniently 
defined observables are consistent with the Standard Model (SM) predictions, thus leaving 
little room for New Physics (NP). Since the non-perturbative QCD (hadronic) uncertainties 
are still large in most cases, the comparison between theory and experiment is not yet at the 
precision level and the quantitative statements about the size of possible NP contributions 
are often subjects of controversies. 

Much of the experimental activity has been devoted to B — > K*£ + £~ decay, for which 
the improvement on theoretical (hadronic) uncertainties is hard to achieve, apart from a 
few asymmetries that will be studied at CMS, LHCb and in the new generation of B- 
physics experiments. In contrast to the decay to K*, a substantial improvement in the 
determination of the hadronic form factors entering the theoretical description of B — > 
K£ + £~ is realistic to expect very soon. In that respect the recently reported result on 
Br(B — > K£ + £~) by BaBar [2] is likely to become a major constraint in the NP searches. 
The ongoing experimental effort to detect another b — > s mediated decay, B s — > , has 
been greatly improved after LHCb was able to set the upper bound on Br(£> s — > n + 
that got very close to the value predicted in the SM [3]. In this paper we discuss how 
these two decay modes can be combined to give us complementary information about the 
potential physics contributions from BSM. 

The most important effect of physics BSM in B s — > / u + fi" is expected to come from 
a coupling to the scalar and/or the pseudoscalar operators. If that scenario is verified in 
Nature, it would strongly affect B — > K£ + £~, whereas the three transverse asymmetries 
4?' im ' re) (g 2 ) of B -> K*£+£~ decay would remain unchanged with respect to their shapes 
predicted at low g 2 's in the SM [6]. 1 If, instead, the NP alters the couplings to the 
semileptonic operators oc sj^Pi^b, then the shapes of the mentioned asymmetries would 
change too. 

In this paper we assume that the NP does not couple with tensor operators (specified 
below). It turns out that this assumption can too be tested by measuring A FB (q 2 ), the 
forward-backward asymmetry in B — > K£ + £~ , but in the region of large q 2, s in which the 
non-zero tensor couplings would entail a large enhancement of A FB (q 2 ). 

In what follows we will show how and why the two decay modes, B s — > and 
B — > K£ + £~, are complementary, and after combining the recent experimental results 
with our current theoretical knowledge of the hadronic matrix elements we will discuss the 
resulting constraints on NP. 

In handling the B — > K£ + £~ decay we computed the form factors by using the simu- 
lations of quenched QCD on the lattice (LQCD), which appear to be compatible with the 
results obtained by evaluating the QCD sum rules near the light cone (LCSR) [8]. These 
results are merely an illustration of the potential of LQCD in pinning down the hadronic 

1 Please note that the asymmetry A^\q 2 ) has been introduced in ref. [7] and since then abundantly 
discussed in the literature. 
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errors in this decay. Very soon these results will be substantially improved by using the 
available QCD configurations that contain the effects of dynamical quarks [9]. Notice that 
the lattice results are more reliable in the larger half of the g 2 -region available from this 
decay (14 GeV 2 < q 2 < 20 GeV 2 ). For that reason it would be desirable to have the ex- 
perimentally established partial decay width of B — > K£ + £~ at g 2 's that are also accessible 
in modern lattice QCD studies. In that way we would be far more confident about the 
quantitative statements made from this kind of analyses. 

The effective Hamiltonian describing the b — > s£ + £~ transitions at low energy is [10] 



i=l i=7,8,9, 10,P,S,T,T5 ^ ' 



where the twice Cabibbo suppressed contributions (oc V ub V* s ) have been neglected. The 
operator basis in which the Wilson coefficients have been computed is [11, 12]: 

7 = ^m b (sa^P R b)F^, 0' 7 = ^m b (sa^P L b)F^ , 

9 9 

Os = -m b (sa^T a P R b)G^ a , 0' s = -m b (sa ^T a P L b)G^ ' a , 

9 = - 2 (s % P L b) (£Y£) , 0' 9 = - 2 {- sl ,p R b) {hn) , 

9 9 

2 2 

10 = e - 2 {s llx P L b){h»ld), G' w = ^(s lfl P R b)(h^ 5 £), 

9 9 

Os = ^P R b) (ft) , Q' s = ^ (sP L b) (ft) , 

2 2 

lD7T z Ion* 
2 2 

O t = -^(sa^b)(£a^£), T5 = -^{sa^b){£a^ l5 £), (2) 

l07T z lbn z 

where Pl,r = (1 =F 7s)/2, £ = e or /x, and the explicit expressions for Ci_6 can be found 
in ref. [11]. Short distance physics effects, encoded in the Wilson coefficients Cj(/x), have 
been computed in the Standard Model (SM) through a perturbative matching of the effec- 
tive with the full theory at fi = mw, and then evolved down to fi — m b by means of the 
QCD renormalization group equations at next-to-next-to-leading logarithmic approxima- 
tion (NNLO) [11]. It is customary to reassemble Wilson coefficients multiplying the same 
hadronic matrix element into effective coefficients appearing in the physical amplitudes, 
namely [13], 

C s + C 3 -lc + 20C 5 -^-C 6 
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where the function Y(q ) is 

4 fi4 64 1 / 4 fi4 

*V) = ^C 3 + -C 5 + -C 6 --h(q 2 ,0)^C 3 + -C 4 + 16C 5 + -a 

+h(q 2 ,m c ) {^d + Ct + dCs + WC^ - l -h(q\m b ) (jC 3 + ^C A + 76C 5 + y tf fl J , (4) 
and 

2 f arctan — = z > 1 

1 m ^ z < 1 

2 



(5) 

with z = 4m 2 ,/ q 2 . To make the notation less heavy, in what follows we will drop the 
superscript "eff" in Wilson coefficients, while tacitly assuming the redefinitions (3). Note 
that in the SM the right-handed flavour violating couplings or to the operators Os,p,t,T5 
are absent and therefore the Wilson coefficients C' 7 8 g 10 = C^ PT ^ T5 = 0. 



2 Physical processes as constraints 

2.1 B s -> £+£- 

One of the most promising decay modes expected to reveal the effects of NP at LHCb is 
B s — > n + ji~ . Using the effective theory (1) the only operator contributing to the amplitude 
of this process in the SM is (9io, and the expression for the branching fraction of this decay 
reads: 



Br (B s -> £ + £-) SM = r Bs ^\V tb V;fm Ba m 2 Mrn 2 B J\C W \ 2 fl , 



G 2 F a 2 



(6) 



where /3i(q 2 ) = a/1 — Amj/q 2 , i = /i or e, and the B s meson decay constant is defined via 

®\^ ll PLb\B.{p))= % -f B j> ll . (7) 

On general grounds the contributions of physics beyond the Standard Model (BSM) can 
modify the above expression to 
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thus lifting the helicity suppression exhibited in the SM expression (6). This is why it 
is interesting and important to experimentally investigate this decay both for i = and 
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Figure 1: Increase of the Br(i? s -> [i+u, ) with 
band defined by la errors on the input parameters. The brightly shaded area depicts the experimentally 
allowed values for the decay mode. C' w — C^) = has been used. 



t = e. In the SM the electron mass severely suppresses the decay rate, and a clean detection 
of the B s — > e + e~ events would be a clear signal of NP effects. In the last formula the 
presence of the right handed (non-SM) couplings would induce C' 10 ^ 0, while the new 
scalar (pseudoscalar) couplings would entail the non-zero Cg (Cp ). 

Several specific NP models that allow for large values of Cs suggested a possibility 
of observing a large number of B s — > fi + fi~ events. Recent experimental activity at the 
LHC [3], however, showed that such a large enhancement does not occur and the current 
upper bound on Br (B s — > fi + fi~) is quite close to the SM value. More specifically, 2 

Br (B s /2 + fi-) cxp < 4.1 x 10" 9 , 

(9) 

Br (B s -> fj+fx-y h ' SM = (3.3 ± 0.3) x KT 9 , 

where for the SM estimate we used the parameters given in the Appendix B of the present 
paper (c.f. tab. 2). Although the desired enhancement by orders of magnitude with respect 
to the SM does not occur, the possibility of seeing the NP signal from this decay mode is 
still alive. One can even envisage the possibility of Br (B s — > smaller than the one 

predicted in the SM, as it can be easily seen from eq. (8). 

If, for the moment, we only consider the possibility that the SM is extended by allowing 
a coupling to the scalar operator, i.e. by keeping C 10 = C^} = 0, then from eq. (8) one 

2 Recently, it has been noted that the effect of B s — B s mixing should be taken into account [4] and 
whose net effect amounts to replacing Br(B s —> /i + /i~) cxp — > (1 — y s )Br(B s —> /i + ^i~) cxp , where y s = 
(AT /2T)b s — (9.0 ± 2.1 ± 0.8) % was also measured by LHCb [5]. This correction is already incorporated 
in the corrected experimental value (9). 
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gets 



\C S ~ C 3 \ < 0.08 (1 a) , (10) 

that we illustrate in fig. 1. This limit is actually not too far away from \Cs — C' s \ M = 0, 
which is why this bound becomes a very severe constraint on the NP models with scalar 
operators. The above bound still allows for an observation of Br(5 s — > e + e~) <2x 1CT 9 , 
unless the lepton flavor universality is not respected by NP. Obviously, the bound (10) 
does not give us any information about Cs and C' s separately. For that we would need a 
complementary information that can be inferred from B — > K£ + £~ . 

Recent research about exploring the constraint from the experimental bound on By(B s — > 
/U + /i _ ) has been reported in ref. [14]. 

2.2 B K£ + £~ 

There are many papers in the literature dealing with this decay. We were able to check 
the expressions given in ref. [15] with which we agree. The full distribution of this decay 
is conveniently expressed as, 

d2T d ^f d ZT = g ^ 2) + h ^ q2) cos 9 + c ^ q2) cos2 ^ (n) 

where for short Te = T(B — > K£ + £~), q 2 = (p e + + Pe-) 2 , and 9 is the angle between the 
directions of B and of l~ in the center of mass frame of the lepton pair. In terms of specific 
Lorentz components, F,i{q 2 ) (i = S, P, V, A, T, T5), the explicit expressions of the functions 
on the right hand side of eq. (11) are: 

aM 2 ) = C(q 2 ) [q 2 (^ 2 )\Fs( q 2 )\ 2 + \F P (q 2 )\ 2 ) + ^ (|i^(g 2 )| 2 + \F v (q 2 )\ 2 ) 

+ Am 2 m 2 B \F A (q 2 )\ 2 + 2m, (m| - m\ + q 2 ) Re {F P {q 2 )F* A {q 2 )) 



b e (q 2 ) = 2 C(q 2 )[q 2 [/3 2 (g 2 )Re (F s (q 2 )F*(q 2 )) + Re (F P (q 2 )F* 5 (q 2 ))] 

+ m t [y/\(f)^(q 2 )Re (F s (q 2 )F^q 2 )) + (m| - m\ + q 2 ) Re (F T5 (q 2 ) F\{q 2 )) ] } , 

dtf) = C(q 2 ) [q 2 W)\F T { q 2 )\ 2 + \F T5 (q 2 )\ 2 ) - {Wa^ 2 )? + \Fy(q 2 )\ 2 ) 

+ 2m £V ^)/3,(g 2 )Re (F T (q 2 )F^(q 2 )) ] , (12) 

where 
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with 



A(g 2 ) =g 4 + m| + m^-2( 
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(14) 
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m B -m K 
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(15) 



In the above expressions we employed the standard decompositions of the hadronic matrix 
elements in terms of the form factors, namely, 



(K(k)\s^b\B(p)) 



(P + k) t 



m B -m~ K 
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ttl b + m K 

After integrating eq. (11) over q 2 = (p — k) 2 , one obtains 

cffVcos 9 



(16) 
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A t + B e cos 9 + Ce cos^ 9 
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where obviously 
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with g^ in = 4m 2 , and g^ ax = {m B — m K ) 2 . For a partial decay width, one would obviously 
choose different q^ , , a Y . Finally, the integration over 6* leads to the full decay width, in 
which the term proportional to cos 9 drops out. The latter survives in the forward-backward 
asymmetry and we have, 



=2(A e + \Ce) , A 
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Figure 2: Br (£? -> if^ - ) = r B T e as a function of \Cs + C' s \. Bright shaded horizontal band corre- 
sponds to the recently measured Br (B — > Kl + £~) at BaBar. We used given in eq. (19), and t b = t b <>, 
to conform with the experimental practice when combining neutral and charged B-modes. 



Clearly these two observables are independent as they involve different pieces of the dis- 
tribution (11). In ref. [15] another useful quantity has been introduced, 



^ 2 {At + Ct) 



F« H = , (20) 

J- e 

which in the SM is proportional to mj, but can receive important contributions in various 
scenarios of NP. Before continuing, we wish to stress that in the SM A e FB = 0, and it 
remains zero even if the Wilson coefficients Cj \ 10 received large NP contributions. Its 
non-zero measurement would be a clean signal of NP and therefore its experimental study 
would be highly welcome. 

For our purpose, it is important to note that all the functions in eqs. (12,15) involve 
the sum of Cj + and therefore the full branching ratio provides us with a constraint 
that is complementary to eq. (10). The recently measured [2] 

Br (B -»■ K£ + £~) exp = (4.7 ± 0.6 ± 0.2) x 10" 7 , (21) 

is compatible with the SM prediction, 

C (7.5 ± 1.4) x 10~ 7 LQCD, 
Br(B^Kl+r) su =\ , (22) 

[ (6.8 ± 1.6) x 10~ 7 LCSR. 

which, at this stage, we take to be 

Br (B -> K£ + £-) SM = (7.0 ± 1.8) x 10" 7 , (23) 
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thus covering all the values allowed by the two methods of computing the form factors, 
LQCD and LCSR. Allowing for non-zero Cg then leads to 

\Cs + C' s \<1.3 (la), (24) 

as illustrated in fig. 2. To obtain eq. (24) we had to assume that NP does not alter 
the SM values of C*}\ which enter the expression for Br (B — > K£ + £~). That assumption 
can be tested experimentally through the study of low-g 2 behavior of three transverse 
asymmetries discussed in ref. [6] which exhibit three very important features: (i) they have 
small hadronic uncertainties, (ii) their shapes are highly sensitive to Cyg, and (iii) they are 

completely insensitive to C^ p . 

Another important comment concerning the constraint (24) is that it is obtained by 
including the la experimental uncertainty, and by using the form factors that are either 
obtained from the numerical simulations of quenched QCD on the lattice (see Appendix B 
of this paper), or in the QCD sum rule analysed near the light cone [8], respectively labelled 
as LQCD and LCSR in eq. (22). LQCD results appear to be consistent with those obtained 
from LCSR. The three relevant form factors, / +j o,t(o 2 ), will soon be improved by the new 
generation of unquenched lattice QCD simulations. Several such studies are underway [9]. 
Notice that the improvement of / +i o,t(o 2 ) is much more realistic to expect than those 
parameterizing the B — > K* transition matrix elements, because the latter decay involves 
many more form factors, including at least three that suffer from very large uncertainties 
(see e.g. ref. [16]). 

3 Constraints on the Scalar (Pseudoscalar) Couplings 

Couplings to the scalar and pseudoscalar operators are particularly interesting in the frame- 
work of supersymmetric (SUSY) extensions of the Standard Model. A very detailed con- 
sideration of the SUSY contributions to B s — > £ + £~ has been made in ref. [17]. The relation 
between the Wilson coefficients Cg P and Cf LLRRLRRl defined in ref. [17], are: 

C s = X (Cf R + C[ L ) , C P = X (C[ R - C[ L ) , 

C' s = X (C RR + C RL ) , C'p = X (C RR - C RL ) , (25) 

where X = ir/{V2G F aV*V tb ). From ref. [1 7] we learn that the SUSY contributions to the 
box and penguin diagrams can modify Cg P as follows: 

• Diagrams with one charged Higgs boson propagating in the box can give non-zero 
contribution via coupling to the left handed parts only. Furthermore, they verify 
Cg + = —Cp + . The right-handed couplings, instead, are suppressed by the strange 
quark mass, C' S H+ = C' P H+ = 0. 

• Diagrams with charginos propagating in the box also give rise to Cg and C R , but 
leave C s x = C' p x = 0. Moreover, if for example the masses of squark and sneutrino 
in the box are degenerate, then one again obtains C$ = —C p . 
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• Z°-penguin diagram, with super-particles propagating in the loop, cannot generate 
a contribution to C^ p due to the vector coupling to Z°. 

• -£f°-penguin, instead, can give a sizable contribution which verifies 

C f = -Cf, Cf = Cf, (26) 

up to the electro-weak symmetry breaking corrections that are proportional to the 
mass splitting between the SUSY (MSSM) Higgs bosons, namely m 2 H0 — m 2 A0 . 

The situation similar to the if °-penguin case above also happens in the models with vector 
leptoquark states. In those models the non-zero contributions to C,s are possible and they 
satisfy the relation Cs = ±Cp, as well as C' s = ^fC'p [18]. On the contrary, the models with 
scalar leptoquarks cannot generate any sizeable value of C^ p . These are obviously only 
two among many scenarios of physics BSM in which the coupling to scalar (pseudoscalar) 
can be generated. 

In the remainder of this section we will combine the constraints discussed in Sec. 2 
and apply them to two particular scenarios: (1) Cf ^ 0, while C[ = C { p> = 0, and (2) 
C^p 7^ 0, while keeping C' l0 = C*g = 0. In each case we will also allow the non-zero relative 
phase between the left- and right-handed couplings (Wilson coefficients). 

3.1 Cf ± 

We first focus on the scenario in which Cs and C' s can be different from zero. Besides 
|CW| ^ o, NP can induce new weak phases, which through our constraints in eqs. (10,24) 
cannot be studied separately. Instead, one can study the impact of the relative phase, 
A(f) S = 4>'s ~ <t>s, by using 

\C S ±C' S \ 2 = \C S \ 2 + \C' S \ 2 ± 2 \C S \ \C' S \ cos(A0 5 ) , (27) 

and explore the possible values of \Cs\ and \C' S \ for various Acps G [0, it] that are compatible 
with the constraints (10) and (24). The bound on B s — > fi + fi~ is much more compelling, 
and for any non-zero relative phase, A05, the constraint provided by B — > K£ + £~ becomes 
essentially superfluous, except in the case of A(p s = when the constraint (10) describes 
a stripe in the plane (|Cs|, \C' S \) that is cut off by the constraint (24). This situation 
is illustrated in fig. 3 for three representative cases, Acps = 0, ir/2, and tc. We see that 
the situation in which A<fis = is indeed peculiar, and only for very small values of the 
relative phase A(f>s the sizable couplings to new physics via scalar operators are possible. 
Otherwise a non-zero relative phase entails a reduction of available \C$ ' |, as we show in 
fig. 4. Two important comments are in order: 

• Any value of \Cs\ and/or \C' S \ allowed by the constraints (10) and (24), for any value 
of the relative phase A<fts, are consistent with the branching fraction of the inclusive 
B — > X s £ + £~ decay, which has been measured at the _B-factories at low g 2 's [19, 20], 
resulting in an average [21], 

r^'^f-^W = 1.6(5) x 10- . (28) 

Jl GeV 2 "5 ex P 
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By using the formulas presented in refs. [21, 22], we obtain in the SM 

f °' V ' dBr <^f-^W = 1.59(17) x 10- , (29) 

Jl GeV 2 HQ th ~ SM 

where, instead of the usual practice to normalize by the inclusive semileptonic T(B — > 
X c ev) decay, we actually use the tree-level decay rate proportional to the fifth power 
of the 6-quark mass, i.e. 

4 a 2 Gl \V tb V*\ 2 mf, , s 7 

B ° = TB 3 ( J)5 " = 3 - 41 ( 47 ) x 10 > ( 3 °) 

which is now possible thanks to the fact that the value of the 6-quark is by now very 
well determined from the multitude of techniques of QCD sum rules and modern 



A0,, = O " Afe = 90 ° A0, = 180 




Figure 3: Constraint on \Cs\ and \C' S \ obtained from Br (B — > Kl + £~) are represented by the brightly- 
shaded area, whereas the one deduced from Br(B s — > is described by the dark shaded region. 
Of course, only the overlap of both regions is consistent with the constraints. It is also consistent with 
Br (B X s £ + £~) and |C^| 7^ it does not modify the transverse asymmetries in B — s- K* ( + £r . Plotted 
are the cases with three specific choices of the relative phase Alps- 
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Figure 4: Maximal value of \Cs\ or \C' S \ allowed by the constraints (10) and (24) as a function of the 
relative phase A(f>s £ [0, 7r]. 
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simulations of QCD on the lattice. Our result agrees very well with the one obtained 
in ref. [21] that was obtained by normalizing to T(B — > X c ev). In practice we take 
the average quark mass quoted by PDG and converted it to the pole mass using the 
NNLO perturbative QCD corrections [23]. 

The scalar contribution to B — > X s fi + fi~ has been computed in ref. [24], and leads 
to a term that adds up to Br (B — > X s £ + £~) su and reads, 

dBr(B -»> X s n + n~) 3B q 2 \ 2 q 2 



4'» 2m 2 



dq 2 

When integrated between 1 and 6 GeV 2 , this leads to a tiny correction, 



I- 6 dBT(B -» Xs» + V-) dq 2 = 1 59(17) x 1Q ~6 ^ + Q 0Q7(1) Q Cg f + j^f)] ^ 
Jl GeV 2 "<? 

(32) 

so that even the largest allowed values for \Cg \, displayed in fig. 4, result in a 
negligibly small correction to the inclusive B — > X s fi + fi~ rate. 

• We stress again that the NP scenario in which only the coupling to the scalar opera- 
tors are allowed would not modify the SM predictions of the g 2 -shapes of three 
transverse asymmetries that can be measured from B — > K*£ + £~ decay, namely 

4 2 ' im ' re V) [6]. 

3.2 C ( p ] ^ 

The situation is slightly more complicated in the case of Cp ± C' P because both phases are 
needed. It is easy to see that eq. (8) gives us 



Up — Up -+- z — ^ — o 10 

m B, 



aG F f Ba m 2 B \V tb V*\ W m Bs T B J^m 2 B ) 



and therefore considering only the relative phase, A(p P = (ppi —<pp, is not enough. Instead, 
we write 

\C P \ 2 + \C' P \ 2 - 2\C p \\C' p \cos(A(Pp) 

+ C 2 + 2 C 10 [\C P \ cos <j)p - \C' P \ cos(0 P + A0 P )] < (0.150) 2 , (34) 

where, for brevity, we use Cio = 2Cf^mi ) m^/m 2 Bs . Besides Acfip we choose to vary the 
phase of Cp. A compact analytical expression similar to eq. (34) that constrains Cp + C' P 
from Br (B — > Kl + l~) cannot be obtained. Instead, we get, 

\Cp + C' P \ mi=0 <1.3, |C P + C$,- 0.33| < 1.3 (la), (35) 
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<pP = °, A<f> P = P = ", A0 P = 90 0j> = Afo = 180 




Figure 5: Allowed values for \Cp\ and \C' P \ correspond to the overlap of the constraints obtained from 
Br [B — > K£ + l~) (light shaded area) and Br (B s — > (dark shaded area). Illustration is provided for 

three various values of the relative phase Atpp and for <fip = 0. 



in the massless and massive lepton case respectively. In what follows we will use rri£ = 
m M ^ 0. 

We first fixed <pp = and varied the relative phase A<pp G [0, tt]. Typical examples 
are shown in fig. 5. We then let (ftp G (0, 27r], and observe that the possible values of \Cp\ 
and \C' P \ consistent with (9) and (21) are smaller than in the (f)p = case. For any fixed 
(f>P the situation is similar to what we observed in the case of \Cg\ and \C' S \ (c.f. fig. 4), 
namely that for larger A0p the possible values of \Cp\ and \C' P \ are smaller than in the 
case A0p = 0. In other words the most space available for NP occurs when the phases of 
\C P \ and \C' P \ are aligned (A0p 0), and even more when the new physics phase 0p ~ 0. 
This is illustrated in fig. 6. 

Furthermore the observations similar to those we made in the end of previous subsection 
apply also in this case: 

• Our result that \Cp\ < 1.0 for any value of 4>p and for any A0p, is consistent with 
the observed branching fraction of the inclusive decay rate (28) which is modified by 
the presence of the pseudoscalar operator in the same way it was in the case of the 
scalars [24], namely, 



dBr(B -> X s ii + ri- 



dq 2 



3B ( g 2 \ 2 g 2 



c ( p } 2ml V ml) ml 



|Cp| 2 + |^| 2 ). (36) 



• Notice also that the non-zero values of C^p ^ cannot modify the SM predictions of 
the low-g 2 shapes of three transverse asymmetries, A^' Mn ' , currently studied in the 
B -»■ K*£+£- decay at LHCb. 

3.3 Peculiar case of C$,p ^ 

Before closing this section, we would like to comment on the case, often discussed in the 
literature, in which the NP can couple via Cs,p ^ but with C' s P = 0. The available 
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range of values for |Cs 5 p| 7^ consistent with the constraints provided by Br (B — > K£ + £ ) 
and Br (B s — > is depicted in fig. 7. As in the previous cases the largest range of 



</> P = 0° 



<p P = 90 " 



ij> P = 180° 




0.1 0.2 0.3 0.4 0.5 

icy 




Figure 7: Allowed values for \Cs\ and \Cp\ obtained by combining the experimental information on 
Br (B — > K£ + £~) (light shaded area) and the upper bound on Br (B s — > (dark shaded area). Illus- 

tration is provided for three various values of the relative phase <pp. 



\Cs,p\ ^ is obtained when the pseudoscalar coupling is real, (ftp = 0. The result is of 
course invariant with respect to the change of the phase <ps ■ A particularly important 
observation that can be made in this case is that the current constraint provided by B — > 
K£ + £~ is redundant, but that situation could radically change if the errors on B — > K 
form factors were significantly reduced. To illustrate that effect we keep the central values 
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of the form factors fixed and reduce the errors by 20%. In that hypothetical situation the 
constraint coming from the measured and theoretically evaluated Br (B — > K£ + £~) are not 
compatible with the SM (within the la accuracy), and therefore B — > K£ + £~ becomes an 
essential constraint to the values of possible IC^pj. The corresponding plots are presented 
in fig. 8, where we only show the cases for which the overlapping region (satisfied by both 
constraints) exists. For cf)p > 40° such a solution would not exist, which would be very 
valuable information about new physics. 



<p r = " 4> P = 20 " 4> P = A0° 




Figure 8: Same as in fig. 7, except that the errors on the hadronic form factors relevant to B — > K£ + t 
are reduced by 20%, and we plot the cases with (ftp = 0, 20° and 40°. 



The above example only further highlights the importance of reducing the errors on the 
B — )■ K transition form factors by using the currently available lattice QCD configurations 
that include iVf = 2, 2 + 1, and even 2 + 1 + 1 dynamical quark flavors. We should also 
stress that the simultaneous experimental effort in measuring the partial decay width of 
B — > K£ + £~ at several moderately large and large values of q 2 would be highly welcome 
because for those momentum transfers the uncertainties of the form factors computed in 
LQCD are under much better control than those at low g 2, s. Effort in that direction made 
in ref. [2] is highly welcome. 

4 Constraints on C^jJ 

In this section we focus on the NP contributions that might arise from the couplings to 
the operator O\o and O' l0 , assuming that the Wilson coefficients C^ P = 0, as in the SM. 3 
Specific realizations of the two-Higgs doublet models have been discussed in great details 
in refs. [17, 25]. Note that our Wilson coefficients CfJ are related to the ones defined in 
ref. [17] as: 

Cio = X (CY R - CYl) , C[ = X (C\[ R - Cl[ L ) , (37) 

where X is the same one defined after eq. (25). From refs. [17, 25] we learn that: 

3 Here we also tacitly assume that the \ T , which enter the expression for Br (B — > K£ + £~), remain 
at their SM values. 
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The Z°-penguin, with a charged Higgs running in the loop, gives rise to 



.2 ^ 

m^tan 2 /?' " m 2 ^ 



C?o° H+ C^V-^tan 2 /3. (3* 



Therefore a non-zero contribution in the scenario is conceivable for either small or 
large tan/3, although the large tan/3 in C' w is suppressed by the strange quark mass. 

• The Z°-penguin, with a gluino running in the loop, is only relevant at larger tan /3 
and the corresponding contributions are such that Cf 9 = C'^ 9 '. 

• Box diagrams are highly suppressed and give no interesting contributions. 

As far as the leptoquark models are concerned, C\q and C' 10 can be non-zero in both classes 
of models, namely with scalar of vector leptoquarks. Interestingly, however, the change in 
C{o implies the change in Cg too. For more details about this issue see ref. [18]. 

In what follows we proceed in a way similar to the previous section, and use eqs. (9) 
and (21) to obtain 

|c 10 + eg < 4.4 , |c 10 - cg < 4.8 (1 a) . (39) 

It is now sufficient to study the impact of the relative phase, A<p = <fi 10 i — </> 10 , by using 

|c 10 ± eg 2 = |cg 2 + \c 10 \ 2 ± 2 |cg |cg CO s(a</>) . (40) 

In fig. 9 we illustrate the resulting constraints for three distinctive cases A<p = 0, 7r/2, and 
7r. We see that for A(f> = the main constraint comes from B — y K£ + £~, whereas in the 
case of A0 = 7r the decisive constraint is B s — y fi + fi~, which is easy to understand from 
eq. (40). In the intermediate case of A</> = ir/2 the two constraints are equivalent. 

We also checked the hypothetical scenario in which the measured Bi(B s — y /i + /U _ ) 
coincides with its value predicted in the SM (9). As a result the allowed region in the 
|Cg-|Cg plane is depicted by the yellow stripe in fig. 9. 

Contrary to the previous section, in this case the inclusive branching fraction (28) 
provides us with a valuable new constraint. The contribution from 0±q to the differential 
decay rate can be extracted from ref. [24], and for the massless lepton pair it reads, 



dBr(B -> X a n+ii 



R / ^2 \ 2 / „2 



dq 2 

and therefore we can write 

GeV~ 



c w m b V m l J V m l 



! + 2-M (|Cio| 2 + |C; o | 2 ) , (41) 



10 6 x/ 6 dBl {B ± f^^ dg 2 = 0.69(9) + 0.058(6) ( |C 10 | 2 + |^ | 2 ) ■ (42) 

J l GeV 2 "9 V J 

When combined with the experimental value (28) we get 

4.7 < Idol 2 + |Cg 2 < 28.9 (1 a) . (43) 
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= 0" Af = 90" Ad> = 180" 




Figure 9: Constraint on \Cxq\ and \C' 10 \ obtained by combining the experimental information on 
Br (B — > Kl + t~) ("light shaded area) and Br(B s —> (c/arJc shaded in the plots). We used the 

hadronic form factors and the decay constant given in Appendix A and B. Three plots correspond to 
three specific choices of the relative phase A<fi = 0iq' — <pio indicated in each plot. The domain inside the 
dashed curve is allowed by the inclusive decay, as indicated in eq. (43). The region satisfying all three 
constraints is within the thick curve. See text for the explanation about the yellow region. 



As before, we account for the la uncertainty around the central experimental value and 
take the lowest/largest possible values to obtain the limits in eq. (43) which describes a 
disc in the |Cio|-|C{ | plane, as shown in fig. 9 (dashed curves). The situation is now more 
interesting as it depends considerably on the value of the relative phase. For A0 = the 
constraint coming from B — > K£ + £~ is overwhelming and the one inferred from B s — > fi + fi~ 
is only marginal. For A0 = it the two constraints exchange the roles, and the most stringent 
constraint comes from B s — > . In the intermediate situation with A0 = 7r/2 the two 
constraints are equivalent and have the shapes similar to that coming from B — > X s n + fi~ . 
In fig. 10 we show the possible values of | C^o I compatible with all three constraints that 
and for any A<p e [0, it]. 

Another difference with respect to the (pseudo-)scalar operators discussed in the previ- 
ous section is that the low q 2 dependence of the three transverse asymmetries, Ay' im ' re ^(g 2 ), 
extracted from the full angular analysis of B — > K*£ + £~ decay, are different from their SM 
shapes when the Cj'J are modified by the NP contributions (see ref. [6] for details). In 



particular the slope of the asymmetry A^ e \q 2 ) is highly sensitive to the value of C\ 



dAf\q 



dq 2 

where R is a convenient ratio of the B — > K* form factors. 4 



'More specifically, R = {V/Ti)/(m B + m K .) « (A 1 /T 2 )/(m B - m K ,) [6]. 
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0' • 

1 2 3 4 5 

ICiol 

Figure 10 : Dark region describing the possible values of \ C\q \ and \C' 10 \ for any relative phase Atfr £ [0, n] , 
obtained from the measured Br (B — > K£ + £~~), the upper bound on Br (B s — > fi + fi~), and the branching 
fraction of the partial inclusive decay rate, c.f. eqs. (28,42). A thick dot corresponds to the Standard Model 
value, Cio = -4.103. 

5 On importance of measuring the forward- backward 
asymmetry in B — > K£ + £~ decay 

The observation that the spectrum of b — > s£ + £~ decays in the region of large q 2 > 
15 GeV 2 's is not plagued by the cc resonances opened numerous possibilities for testing 
theory against experiment [26]. This is the region in which a major progress in taming 
the hadronic uncertainties by means of LQCD is possible and therefore a more reliable 
extraction of physics BSM from experiment should be possible. 

One quantity that we find particularly interesting to study at large g 2 's is the forward- 
backward asymmetry Ap B (q 2 ). In the SM this quantity is zero and remains as such even 
if the NP considerably modifies the values of the Wilson coefficients 10 . On the other 
hand, if the NP gives rise to the new (non-SM) Dirac structures, Ap B (q 2 ) can quite ap- 
preciably differ from zero. To our knowledge this observation was made for the first time 
in ref. [27]. The expression for Ap B (q 2 ) used in ref. [27], however, differs from the one 
reported in ref. [15]. We checked both formulas and agree with the one given in ref. [15]. 

To get a better insight in the impact of NP on A £ FB (q 2 ), we expand it in powers of the 
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lepton mass and write 5 



C(q 2 ) m B - m K 



y/Mcf) fo(q 2 ) <2(C S C T + C P C T5 ) q 2 f T (q 2 



+ me 



C s C 9 (m B + m K )f + (q 2 ) + 2m b (C S C 7 + 2C T5 C W ) f T {q 2 ) + 0(m 



(45) 



All the quantities in the above formula have already been defined in Sec. 2.2. Measuring 
A e FB (q 2 ) at large g 2 's would be highly beneficiary for our quest for NP at low energies. A 
separate experimental study of A e FB (q 2 ) and A FB (q 2 ) would help us discern the first term 
from the second in (45). Notice that the first term is non-zero only if the NP coupling to 
a tensor operator is allowed. Therefore this quantity can be used to test the assumption 
we made in the previous sections of this paper when discussing Bi(B — > K£ + £~), namely 
that the Ct,t5 = 0, as in the SM. 

Moreover, from the inclusive branching fraction (28) to which the tensor operators 
contribute as [24], 



8£? 



HIT 



nit 



\Ct\ 2 + |Ct5| 2 ) 



(46) 



one gets, 



6 GeV dBi(B ->■ X„n + n j , 



1 GeV 2 



dq 2 



dq 2 = 1.59(17) x 10" 6 [1 + 0.66(9) (|C T | 2 + |C7 T5 | 2 )] . (47) 



Contrary to the (pseudo-) scalar case in which the factor multiplying new Wilson coeffi- 
cients is very small, c.f. eq. (32), the corresponding factor multiplying the tensor Wilson 
coefficients is much larger and consequently the constraint provided by eq. (28) is much 
stronger, 



|CV| 2 + \C T5 \ 2 < 2.6. 



(48) 



For the sake of illustration we take CV = CV5 = 1.1 and plot A FB (q 2 ) in fig. 11 at large 
g 2 's by taking C$p = 0. With that choice, and due to the fact that we consider the 



5 For simplicity, here we take the Wilson coefficients to be real. If they were all fully complex then the 
expansion in eq. (45) would look as: 

A FB (q 2 ) = C{ f )mB ~ mK VW)h{q 2 ) Rc{2 [(Cs + C' S )C* T + (Cp + C" P )C* T5 ] q 2 f T (q 2 ) 
+ m e [((C s + C' s )(C g + (%)*) (m B + m K )f + (q 2 ) + 2m b ((C s + C' S )(C 7 + C 7 y)f T (q 2 ) 
+ 2 ((C w + C' W )C* T5 ) f T (q 2 )] +0(m})\. 
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Figure 11: Forward-backward asymmetry in B —> Ku, + fi~ decay. 
Ct — Ct5 = 1-1 [consistent with eq. (48)] and by keeping Cs,p 
obtained with the same Ct — Cts = 1.1, but with Cs = 0, Cp 
Cs — (thin dashed curve). 



A 
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The full curve is obtained with 
0. The dashed curves, instead, are 
1 (thick dashed curve) or Cp = 1, 



decay to the pair of muons, the A FB (q 2 ) ^ everywhere and is strongly enhanced near 
<?max = ( m B -%) 2 ' We then switched Cs or Cp, to illustrate the case when the first term 
in eq. (45) is non-zero. If indeed realized in nature, this latter situation would be relatively 
easy to check experimentally. 



6 Summary 

In this paper we showed that B s — > n + n~ and B — > K£ + £~, the two actively studied 
decays in the 5-physics experiments, provide us with complementary information about 
the potential NP contributions. While the decay amplitude for B s — > is proportional 

to the difference between the Wilson coefficients of the operators of opposite chirality, the 
B — > K£ + £~ involves the sum of these Wilson coefficients. 

We checked the situations in which the NP enters either via the scalar, pseudoscalar 
or the semileptonic operators. To decide which situation is verified in Nature (if any), 
useful information can be obtained from the low-g 2 shapes of the transverse asymmetries 
in B — > K*£ + £~ decay. Those asymmetries are being studied in experiments and have 
an important advantage that the relevant hadronic uncertainties are small. A non-zero 
coupling to the scalar and/or pseudoscalar operator would not modify the low-g 2 shapes 
of these asymmetries. From B s — > fi + fi~ and B — > K£ + £~ we find the absolute bounds, 

l4 /} |<0.7, |C^| < 1.0, (49) 

that are valid for any value of the NP phases. In fact, the values for |C^p| can get 
considerably reduced if the non-zero NP phases are allowed. 
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In the case the coupling to the semileptonic operators O' 10 is modified by the presence 
of NP particles, the transverse asymmetries in B — > K*£ + £~ would have peculiar shapes, 
different from those predicted in the SM. From our study of B s — > and B — > K£ + £~ 

we obtain that 

2-2 < v / |Cio| 2 + |Q P < 4.8, (50) 

regardless of the value of the relative NP phase. Note that the lower bound is fixed by the 
experimentally measured partial decay rate of the inclusive B — > X s fi + fi~ decay. 

In considering B — > K£ + £~ we ignored the contributions from the tensor operators. 
That assumption can also be experimentally tested by measuring the non-zero forward- 
backward asymmetry in B — > K£ + £~ decay. 

Our approach of considering a pair of Wilson coefficients at a time is orthogonal to 
the global fit approach adopted in many recent works. We have checked explicitly our 
results with the results of Bobeth et al in [14] where a fit to real Wilson coefficients CV^io 
included also the experimental observables related to B — > K*£ + £~ and B — > K*'-/. Our 
results in the scenario with complex C\$ and C' l0 (fig. 10) agree well with their presented 
range of real Ciq. We have also checked our allowed regions for Wilson coefficients against 
the results of the global fit presented in Altmannshofer et al in [14]. 

In this paper we focused on the quantities that either have small hadronic uncertain- 
ties or those for which the hadronic uncertainties are likely to be improved soon. This is 
particularly the case with the B — )■ K form factors, that we computed in the quenched ap- 
proximation of QCD, which will soon be improved by including the effect of light dynamical 
quarks. A detailed experimental information about the partial decay rate of B — > K£ + £~ 
in the upper range of g 2 's will become particularly useful because the results for the form 
factors computed on the lattice at larger q 2, s are more reliable and have smaller errors. 



Note Added: 

While this paper was in writing the new results for Br(B — > K£ + £~), measured at LHCb, 
appeared in ref. [38]. Their value is lower than the one reported by BaBar [2], and the 
agreement with the current theoretical estimate of the same quantity is only at the 2<r-level. 
Once the more reliable estimate of the B — >• K form factors computed on the lattice become 
available, we will repeat the analysis presented here by including the 1, 2, 3cr effects. 
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Appendix A: B — > K Form Factors 



The results for the form factors used in this paper are obtained from the analysis following 
the same procedure as the one explained in detail in ref. [28], but by using the (quenched) 
gauge field configurations obtained at finer lattice spacing [a -1 = 3.8(1) GeV]. Those 
configurations have been used to compute the B — > K*j form factors in ref. [29] and we 
refer the reader to that paper for lattice details. 

Besides the form factors f+(q 2 ) and fo(q 2 ) computed along the lines explained in 
ref. [28], we also computed the tensor form factor /r(<? 2 ) appearing in eq. (16). Note 
that the tensor density depends on the renormalization scale that we have set to fi = nib, 
the same scale at which the corresponding Wilson coefficients have been computed. 

It is easy to extend the parameterization of ref. [30] to include the form factor /t(<? 2 ) 
and keep the minimal number of parameters needed to describe the q 2 dependence of all 
three form factors. In terms of poles exchanged in the t-channel, the g 2 -dependence of 
/r(g 2 ; /i) is driven by the states with J p = 1~. 6 The lowest such a state is B*, whose 
couplings to the vector and tensor bilinear quark operators are defined via 

(0\s^b\B*(p,6 r ))=e*; mB *f^, 

(0\5(T lu ,b\B'(p,e r )) = i (p,e*J- Pu e*;) fi.fr) . (51) 
The nearest pole contribution then reads, 

Q m B* 

I/,,.,,, m 2 -m 2 \ m B* /b*9b*bk 



q 2 — m 2 B * 



po lc(2) = -\^fl t 9B t BK 

q 2 — m B » 

where we used the standard definition (B*(e r )\BK) = 9b* bk (k ■ s r ). Similarly, for the 
matrix element of the tensor quark operator we have 

.pole, 2 \ -\(jn B + m K )fl t (u)g B *BK 

=>• /t (? ; n = - 2 r ' 53 

q l — m B » 

and therefore the residua of these two form factors are related, i.e. 

Res f T (q ;n) = -y— Res f+(q ) . (54) 

1 2 -* m B t m B* JB* 1 2 -*™ 2 b* 

It is quite remarkable to note that the pole dominance, that is expected to be a reasonable 

6 Notice in particular that couplings to 1 + states are ruled out when both external states are pseu- 
doscalars. 
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Figure 12: Comparison of the q 2 -dependence of the B — > K transition form factors as obtained in 
quenched lattice QCD with the predictions based on using the QCD sum rules near the light cone. 



at large g 2 's (small recoils), leads to a simple proportionality relation between /t(<? 2 ) and 
f+(q 2 ), namely 



f T (q 2 ; A*) rn B + m K fE* (p) 



f + (q 2 ) m B * f%, 

which is essentially verified in the low q 2 region (large recoils), as obtained in ref. [31]: 

f T {q 2 ) m B + m K 



(55) 



f+(q 2 ) m B 

On the basis of eq. (54) it is judicious to define the scale independent 



(56) 



fV 

m B* JB* 

m B + m K (fi) 



f T (q 2 ) = ^T^HM^/i), (57) 
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which then can be easily included in the parameterization of ref. [30] as 

cfl — a) „ . c(l — a) 



fo(q 2 



l- q y((3m%] 
h{q 2 ) = 77" 



c(l — OLt) 



q 2 /m 2 B *)(l — aq 2 /m B *) 



q 2 /m 2 B *)(l - a T q 2 /m 2 Bt ] 



(58) 



so that only one new parameter («t) is needed to fit all three form factors. The values 
of form factors extracted at various q 2 from our lattice computation are listed in tab. 1, 
where we also used fB*( m b)/fB* = 0.91(3), computed on the same lattices. In terms of the 
above parameters we have 



/+(0) = / (0) = 0.33(4), / T (0) = 0.31(4) 
a = 0.72(14), a T = 0.67(15), = 1.35(15) 



(59) 



which are obtained by either extrapolating the parameters to the I?-meson mass (M-l), or 
by fitting the results from tab. 1 (M-2) to the parameterization (58). 

We should stress that the results for the form factors presented and used in this paper 
are obtained in the quenched approximation of QCD and that they will be updated very 
soon by using the available gauge field configurations in which the effects of the light sea 
quarks has been included. We also note that the results for the form factors presented 
here are compatible with those obtained by using the LCSR which are parameterized as 
follows [81: 



fo(q 2 



0.331(4) 
-g 2 /6.12 2 ' 



U(q 2 ) 



fAq 2 ) 

and illustrated in fig. 12. 



0.162(21) 
1 - g 2 /5.41 2 

0.198(25) 

l-g 2 /5.41 2 ' (l- g 2/ 5 . 41 : 



0.161(21} 



+ 



+ 



\2 ' 



0.173(22) 
4 - g 2 /5.41 2 



,2 ' 



(60) 



q 2 [GeV 2 ] 


fo(q 2 ) 


U(q 2 ) 


f T (q 2 ) 


13.0 


0.51(6) 


0.90(9) 


0.80(11) 


14.5 


0.54(6) 


1.04(9) 


0.92(12) 


15.9 


0.57(7) 


1.21(11) 


1.07(14) 


17.4 


0.61(7) 


1.45(13) 


1.27(17) 


18.8 


0.65(7) 


1.75(18) 


1.52(20) 


20.3 


0.70(8) 


2.19(22) 


1.89(21) 



Table 1: B — > K£ + £ form factors at several values of q 2 . The renormalization scale for the 
tensor form factor is \i = m&. 
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Appendix B: Numerical values of the quantities used 
in this work 

The values of all quantities used in this work are listed in table 2. Two comments are in 
order. 

• We use tb = t b o to respect the experimental practice when combining the charged 
and neutral Bi(B — > K£ + £~) decay modes. 

• Taking the continuum results of three unquenched simulations from ref. [33] in the 
quadrature, one gets the average f Ba = 0.234(6) GeV, which is marginally compatible 
with a spectacularly accurate result reported in ref. [34]. We inflated the error to 
include the central value of ref. [34]. 





-0.304 


[11, 12] 




0.10566 GeV 


[35] 


Cf(m b ) 


4.211 


[11, 12] 


m h (m h ) 


4.25(12) GeV 


[35] 




-4.103 


[11, 12] 


m K 


0.495 GeV 


[35] 


T B a 


1.47 ps 


[35] 


m B 


5.279 GeV 


[35] 


T B o 


1.52(1) ps 


[35] 


m Bs 


5.366(1) GeV 


[35] 


a = a cra (M 2 z ) 


1/128.94(3) 


[36] 




0.234(10) GeV 


[33, 34] 


G F 


1.1664 x 10- 5 GeV -2 


[35] 


\v tb v t * s \ 


0.0403(8) 


[37] 



Table 2: All the scale dependent quantities are obtained in the MS renormalization scheme. The Wilson 
coefficients are evaluated at the next-to-next-to leading logarithms. a s (M z ) = 0.118. 
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